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Abstract: 

Binary systems, in particular eclipsing binaries, are essential sources of our knowledge of the funda- 
mental properties of stars. The ages of binaries, members of open clusters, are constrained by their 
own fundamental properties and by those of the hosting cluster. The ages of eleven open clusters 
are here found by constructing models for the components of twelve eclipsing binaries. The difference 
between the ages we found and the ages of the clusters derived from isochrone fitting is up to 40 %. 
For the binary system V497 Cep in NGC 7160, the difference is about 100 %. Binary systems whose 
primary component is about to complete main-sequence life time, such as V453 Cyg and V906 Sco, 
are the most suitable systems for age determination. Using model results for these stars, we derive an 
expression for sensitive and uncomplicated relative age determination of binary systems (age divided 
by the main-sequence life time of the primary star). The expression is given as logarithm of radii 
ratio divided by a logarithm of mass ratio. Two advantages of this expression are that (i) it is nearly 
independent of assumed chemical composition of the models because of the appearance of the ratio of 
radii, and (ii) the ratios of radii and masses are observationally much more precise than their absolute 
values. We also derive another expression using luminosities rather than radii and compare results. 
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1 Introduction 

Nature is full of signals displaying its past and present 
structure. Astrophysics is extremely specialized in de- 
tecting and extracting such signals from very noisy 
backgrounds. In this context, research on deep regions 
of stars is a real adventure. Our knowledge of the in- 
ternal structure and evolution of stars is mostly based 
on model computations, which take into account and 
combine outcomes of almost all branches of modern 
physics. To test how well our stellar models repre- 
sent real stars, very accurate observational constraints, 
such as mass, luminosity, radius, chemical composi- 
tion, and rotational properties, are required. These 
constraints are obtained from observation of the outer 
regions, which depend on the physical conditions in 
the deep interior. 

Time is one of the essential components of our uni- 
verse. Understanding the evolution of the far and near 
Universe depends on how successful we are in deter- 
mination of the ages of astrophysical objects. In this 
respect, every type of age determination is crucial for 
completing the picture. In this context, we consider 
the observed detached eclipsing binaries (DEBs) that 
are members of open clusters, to obtain some basic 
information pertaining to stars. 

In the literature, thirteen well observed DEBs in 
twelve open clusters are found. Among these sys- 
tems, V818 Tan in Hyades has already been investi- 
gated in detail by the author and his colleagues (Yildiz 
et al. 2006) and is therefore is not included here. 
Thus, twelve DEBs in eleven open clusters are stud- 
ied, namely, V453 Cyg in NGC 6871, V1229 Tan (HD 
23642) in Pleiades, V578 Mon in NGC 2244, V906 Sco 
in NGC 6475, V497 Cep in NGC 7160, V381 Car (HD 



92024) in NGC 3293, V392 Car in NGC 2516, V1034 
Sco in NGC 6231, DW Car in Cr 228, GV Car in NGC 
3532, V615 Per and V618 Per in NGC 869. The masses 
of the binaries vary between 1.54 Mq and 16.84 Mq; 
therefore, this study essentially deals with the struc- 
ture and evolution of the early-type stars. 

The key point for stellar evolution is the increase 
in mean molecular weight, which causes contraction 
in nuclear core. The temperature of the core is in- 
creased by the heat released as a result of contraction, 
and hence nuclear reactions are intensified. This is the 
basic mechanism that controls the evolution and deter- 
mines the time dependence of stellar structure. By 
constructing rotating (for the rapidly rotating compo- 
nents of DW Car) and non-rotating evolutionary mod- 
els for the component stars in the DEBs, we obtain in- 
formation about the age and chemical composition of 
the binaries and hence about their clusters. The infor- 
mation about age and chemical composition is deduced 
by fitting accurate values of luminosities and radii of 
component stars, if possible. We also obtain some sim- 
ple expressions based on the slope of mass-radius and 
mass-luminosity relations for the binaries in order to 
find the age of binary systems easily. 

The ages of stars are derived by means of various 
astrophysical techniques including (i) observed surface 
properties (see, for example, Kharchenko et al., 2005; 
Yildiz et al., 2006); (ii) seismic properties (in partic- 
ular for the stars with the solar-like oscillations, see 
for example, Eggenberger et al., 2008; Yildiz 2007, 
2008); (iii) rotational properties (Mamajek & Hillen- 
brand, 2008); and (iv) radioactive elements (del Peloso 
et al., 2006). The most widely used technique is based 
on the observed surface properties, which is also the 
method of the present study. 
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Figure 1: Observational HR diagram for the DEBs 
in open clusters. 

The Sun, with its precise seismic and non-seismic 
constraints, is the first target to model for an upgraded 
code describing stellar interiors. The agreement be- 
tween the solar models and observed constrains gives 
an idea of the quality of the code. One can ob- 
tain the initial hydrogen (Xo) and heavy element (Zo) 
abundances and the mixing-length parameter of the 
Sun from the calibration of solar models, and use them 
when these quantities are not constrained for any star. 
In Yildiz (2008), Xo, Zo and a are found as 0.70975, 
0.016 and 1.82, respectively, for the solar model with 
heavy element mixture of Asplund, Grevesse & Sauval 
(2005) and with an age of 4.6 Gyr. The details of the 
code used in the model computations can also be found 
in Yildiz & Kiziloglu (1997) and Yildiz (2007, 2008). 

The remainder of this paper is organized as follows. 
In Section 2, we present the observational fundamen- 
tal properties of the DEBs in open clusters available 
in the literature. The models of their components are 
given in Section 3, and the methods of age determi- 
nation based on the slopes of mass-radius and mass- 
luminosity relations are presented and discussed in 
Section 4. Error analysis for age is presented in Sec- 
tion 5. Finally, we list some concluding remarks in 
Section 6. 

2 Observational Data of the 
Binaries in Open Clusters 

Mass transfer obscures the dependence of the funda- 
mental properties of component stars on secondary 
causes, such as, time, chemical composition, rotation, 
and magnetic field. Therefore, semi-detached and con- 
tact binaries are beyond the scope of the present study. 
Accordingly, DEBs, member of open clusters, are stud- 
ied. We have found in the literature studies of twelve 



DEBs in eleven open clusters. The observed proper- 
ties of their component stars are listed in Table 1. The 
sixth and the seventh columns list the ages of the clus- 
ters given in WEBDA database 

(http://obswww.unige.ch) and Kharchenko et al. (2005), 
respectively. For some of the systems, L = iivR^aT^ff 
is not satisfied (see Section 3 and Table 3). Fig. 1 
shows the observed Hertzsprung- Russell diagram (HRD) 
for the component stars of these eclipsing binaries. All 
the stars are either main-sequence (MS) stars or very 
close to the MS. The mass interval of the stars ranges 
from 1.54 (GV Car B) to 16.84 Mq (V1034 Sco A) 
and it covers a large part of the MS in the HRD from 
O-type to F-type stars. The ages of clusters vary be- 
tween 5 My and 310 My. In the upper part of the 
MS in Fig. 1, from left to right, the ages of stars 
increase as theoretically expected. This result shows 
that the time dependence of effective temperature is 
more significant than any other second order depen- 
dence (e.g., chemical composition, rotation), at least 
for these stars. 

The eleven clusters hosting the studied eclipsing 
binaries are at distances ranging from 150 (Pleiades) 
to 2327 pc (NGC3293). Number of stars in the clus- 
ters observed by UBV photometry is good enough to 
estimate the ages of the clusters by isochrone fitting. 
The WEBDA ages taken from different studies vary. 
Kharchenko et al. (2005) essentially use the isochrones 
of Girardi et al. (2002) for age determination. 

These binaries have attracted attention in recent 
years. None of them is among the well known binaries 
in Andersen (1991). However, six of them (V453 Cyg, 
V578 Mon, V906 Sco, V392 Car, V1034 Sco and DW 
Car) are listed in the recent compilation by Torres et 
al. (2010). Here we shall compare the basic properties 
of the stars in Table 1 with that of given in Andersen 
(1991). 

The slope of the mass-luminosity relation derived 
using the data in Table 1 is 3.89. This means that 

LocM^■'*^ (1) 

For the same mass range, we find L oc M^"^ for the 
DEBs with the most precise data on the fundamental 
properties (Andersen 1991). These two values for the 
power of M are very close to each other. The slope of 
the mass-radius relation of the stars in Table 1 is found 
as 0.704. That is to say R oc M° ™. For the well- 
known eclipsing binaries (Andersen 1991), we find that 
R oc M°-^^. The slopes we find for the mass-radius 
and the mass-luminosity relations for binaries given 
in Table 1 are slightly less than the slopes found from 
the well known binaries. The reason may be that many 
stars in the log{A'I /Mq) interval [0.1,0.4] in Andersen 
(1991) have larger radii than their terminal-age MS 
(TAMS ) radius. If we take the mass range as [2.5 Mq, 
16.84 Mq] rather than [1.54 Mq, 16.84 Mq], the slopes 
of the mass-luminosity and the mass-radius relations 
for the well known binaries are found as 3.94 and 
0.723, respectively. These values are very close to the 
values we find for the eclipsing binaries in the open 
clusters. 
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Table 1: The fundamental properties of the double- lined eclipsing binaries, which are members of open 

clusters. 

Star M/Mq R/Rq \og{Te{K)) \og{L/LRQ) twjy) tKh{y) Cluster / 
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.041 
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14.36 ±0.20 
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0. 
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3, 
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V906 Sco A 
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4, 
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.400 X 10*^ 


1, 
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.144 
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X 
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5.39 ±0.40 


2.92 ± 
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4, 
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.860 
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15.0 ±3.0 
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.080 
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0, 
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1, 
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16.84 ±0.48 


7.45 ± 


0, 


.07 


4, 


.471 


± 


0. 


.015 


4. 


.580 
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0, 


.200 


6, 
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6, 


.46 


X 


10^ 
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9.38 ± 0.27 


4.18 ± 


0. 
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4. 
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0, 
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3, 
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0, 
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± 


0, 
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X 
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,875 


± 


0, 


.300 

















1) Southworth et al. (2004a), 2) Southworth et al. (2004b), 3) Southworth et al. (2005), 4) Hensberge et 
al. (2000), 5) Alencar et al. (1997), 6) Yakut et al. (2003), 7) Preyhammer et al. (2005), 8) Debernardi 
& North (2001), 9) Bouzid et al. (2005), 10) Southworth feClausen (2007), 11) Southworth feClausen 
(2006) 



") The values given for the luminosities, radii and effective temperatures do not obey the luminosity 
equation L = AnR'^aT^ (see Section 3 and Table 3 ). 
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3 Ages of the Clusters from 
Models of the Component 
Stars in the Eclipsing Bina- 
ries 

Thw age of a binary system can be found from the so- 
lution of two equations based on the coevahty of the 
component stars with two unknowns. One of the equa- 
tions is written for the difference between the ages of 
the two components of the binary based on observed 
luminosities: 

AtL = tLA - tLB, (2) 

where Jla and tLB arc ages found from luminosities of 
primary and secondary components, respectively. The 
other equation basewd on the radii: 

Afa = tuA — tjiB, (3) 

where Ira and ^rb are ages from the radii of primary 
and secondary components, respectively. Prom the eis- 
sumption of cocvality, AtL ~ Ata — 0. Using nu- 
merical derivatives of A^l and Atn. with X and Z, we 
solve these equations and find X and Z. In general, we 
find that difi^erent combinations of X and Z for a given 
binary give similar results; thus, the solution is not 
unique. If it is not possible to obtain an acceptable 
solution to Eqs. (2) and (3) for a binary system, then 
we either use one of these equations or fit the model 
of one of its components to observed values. 

In the following subsections, interior models of the 
component stars in Table 1 constructed using the stel- 
lar evolution code described in Yildiz & Kiziloglu (1997) 
and Yildiz (2007, 2008) are presented. Prom these 
models, the ages and chemical composition of the bi- 
naries are obtained if possible. 

3.1 On the Fundamental Properties 
of Components of V615 Per and 
V618 Per 

In this subsection, we consider the binary systems V615 
Per and V618 Per in more detail. Their components 

give very different ages. The reason for this differ- 
ence shoul be investigated. For V615 Per, we com- 
pare its components with those of V618 Per and find 
that the luminosity of V615 Per B is lower than is 
exppected for its mass. If we assume that the power 
of mass is 4 for the slope of the mass-luminosity re- 
lation, then log(Ls/L0) = 2.02, which is larger than 
the value found from observations with by amount of 
0.2. Since the slope of the mass-radius relation for 
the components of V615 Per and V618 Per seems rear- 
sonable, we deduce that the effective temperature of 
V615 Per B should be larger than the value found by 
Southworth, Maxted & Smalley (2004a). As a matter 
of fact, the results they give for the effective tempera- 
tures of V615 Per B {Mb = 3.18±0.05 Mq) and V618 
Per A (Ma = 2.33 ± 0.03 Mq) are the same, despite 
the very different masses. 

Reasonable models for the components of V618 Per 
are obtained with Z=0.01 and X=0.735. In the HRD, 



the positions of the models at the ages of the cluster 
given in WEBDA and by Kharchenko et al. are the 
same for V618 Per A. For V618 Per B, the age=19.1 
My (Kharchenko et al. 2005) is more suitable than 
11.7 My (WEBDA). At 11.7 My, V618 Per B is in the 
pre-MS stage. The ages and the chemical compositions 
we find are listed in Table 2. 

The accurate values of the components of V615 
and V618 Per given in Southworth (2004a) do not 
exactly obey L = A-rrR^aT^g. According to South- 
worth (private communication), this arises because the 
same distance for the two components of a DEB is 
enforced in the light-curve analysis. For the system 
V615 Per, however, the problem is with V615 Per 
B, particularly with its effective temperature. For its 
log(L/Lc5) = 1.82 and R/Rq = 1.903 values, its effec- 
tive temperature should be 12000 K rather than 11000 
K. The suggested values are given in bold style in Tar- 
ble 3. 

For age determination of these systems, we use ob- 
servational radii of the component stars as constraints 
for their interior models, because of the complicated 
luminosities and effective temperatures. The age of 
these systems is roughly estimated as 20 My. V615 
Per B seems to be a MS star at this age. 

3.2 On the Fundamental Properties 
of Components of V906 Sco 

The difference between the masses of the components 
of V906 Sco is not very large but the primary compo- 
nent is a TAMS star and therefore very appropriate 
for age determination. The ages are found from the 
solution of Eqs. (2) and (3), based on the coevality 
of the component stars. We solve these equations and 
find X and Z. Using different reference models, we ob- 
tain very different combinations of X and Z, hence the 
solution is not unique. However, the ages found from 
luminosities and radii of V906 Sco A and B are very 
close to each other. For the three cases given in Table 
4, the age found from _Rb (the last column) is slightly 
less than the others. The mean value of ages is 252 
My. Without Irb it is 256 My. 

3.3 V578 Mon 

This system has the lowest values for slopes of the 
mass-luminosity and the mass-radius relations among 
the systems with the early-type components: 

Im = log(LA/LB)/log(MA/MB) = 2.6, 
rm = log(i?A/i?B)/log(MA/MB) = 0.55. 

This means that either the observational values of 
La and Ra are very low or the the observational values 
of Lb and Rb are very high, compared to the values 
expected for their masses. The influence of this point 
on the age determination of the system is crucial. 

Abundances of some heavy elements (including C, 
N and O) are found from disentangled spectra by Pavlovski 
and Hensberge (2005). The oxygen is underabundant 
by about -0.10 dex. If we talce the solar Zs = 0.0122 
(surface heavy element abundance), then Zo of NGC 
2244 is about 0.01. 
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Table 2: The chemical composition and age of the binary systems derived from the calibration of models. 
The ages we find are listed in the fourth column. For comparison, the ages of the clusters given in the 
WEBDA database and by Kharchenko et al. (2005) are also presented. The seventh column reports the 
MS life time of the primary components of the binary systems. In the eighth column, the age found from 
the ratio of radii of the component stars is shown (see Section 4) . The ages found from the fitting formula 
given in Eqs. 6-9 are listed in ninth column 
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.172 X 10 ^ 


1. 


.91 


X 


10^ 


1 
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.91 
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V453 Cyg 
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.7046 


0, 


.0117 
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.05 X 10^ 


9, 
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9, 


.77 
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10^ 


1, 


.18 
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10^ 


9 


.76 


X 


10^ 


9, 


.60 


X 


10^ 


NGC 68' 


V1229 Tau 


0, 


.6946 


0, 


.0158 


1 


.93 X 10^ 


1, 


.5 X 10^ 


1, 


.35 


X 


10^ 


6 


.88 


X 


108 


1 


.24 


X 


108 


1, 


.44 


X 


108 


Pleiades 


V578 Mon 


0, 


.7098 


0, 


.0160 


3 


.50 X lO'^ 


7, 


.870 X lO'^ 


5, 


.01 


X 


10*^ 


1 


.29 


X 


10^ 


















NGC 22^ 


V906 Sco 


0, 


.7200 


0, 


.0170 


2 


.52 X 10^ 


2, 


.400 X 10^ 


1, 


.66 


X 


108 


2 


.53 


X 


108 


2 


.53 


X 


108 


2, 


.53 


X 


108 


NGC 64: 


V497 Cep 


0. 


.7440 


0, 


.0125 


2 


.30 X 10^ 


1, 


.896 X 10^ 


4, 


.57 


X 


10^ 


4, 


.95 


X 


10^ 


2 


.17 


X 


10^ 


2, 


.03 


X 


10^ 


NGC 71( 


V381 Car 


0, 


.7380 


0, 


.0258 


1 


.05 X 10^ 


1, 


.032 X 10^ 


8, 


.71 


X 


10*^ 


1 


.37 


X 


10^ 


1 


.05 


X 


10^ 


0, 


.43 


X 


10^ 


NGC 32! 


V392 Car 


0, 


.6837 


0, 


.0178 


1 


.23 X 10« 


1, 


.123 X 10* 


1, 


.20 


X 


108 


1 


.04 


X 


109 


1 


.41 


X 


108 


3 


.01 


X 


108 


NGC 25: 


V1034 Sco 


0, 


.6961 


0, 


.0226 


5 


.51 X 10^ 


6, 


.967 X 10^ 


6, 


.46 


X 


10*^ 


1 


.02 


X 




5, 


.36 


X 




4, 


.21 


X 


10^ 


NGC 62; 


DW Car 


0, 


.6934 


0, 


.0106 


5 


.22 X 10** 


6, 


.761 X 10^ 


4, 


.79 


X 


10^ 


1 


.56 


X 




4, 


.88 


X 


10^ 


5, 


.72 


X 


10** 


Cr 228 


" (Rotat.) 


0, 


.6951 


0, 


.0078 


5 


.29 X 10^ 


6, 


.761 X 10^ 


4, 


.79 


X 


10*5 


1 


.52 


X 


10^ 


4 


.75 


X 


10^ 


5, 


.75 


X 


10^ 


Cr 228 


GV Car 


0, 


.7046 


0, 


.0104 


3 


.22 X 10*^ 


3, 


.105 X 10^ 


2, 


.82 


X 


108 


4, 


.22 


X 


108 


3 


.21 


X 


108 


3, 


.04 


X 


108 


NGC 35; 



Table 3: The modifications in the fundamental properties of the double lined eclipsing binaries due to 
inconsistent results given in the literature are presented in bold style 



Star M/Mq R/Rq logT,s{K) log(L/La) 



V615 Per A 


4.075 


± 


0.055 


2.291 


± 


0.141 


4.176 


± 


0.033 


2.370 


± 


0.080 


V615 Per B 


3.179 


± 


0.051 


1.903 


± 


0.094 


4.079 


± 


0.045 


1.820 


± 


0.100 


V618 Per A 


2.332 


± 


0.031 


1.636 


± 


0.069 


4.041 


± 


0.033 


1.546 


± 


0.140 


V618 Per B 


1.558 


± 


0.025 


1.318 


± 


0.069 


3.903 


± 


0.045 


0.805 


± 


0.080 


V453 Cyg A 


14.36 




:0.20 


8.551 


± 


0.055 


4.424 


± 


0.019 


4.513 


± 


0.210 


V453 Cyg B 


11.11 


± 


:0.13 


5.489 


± 


0.063 


4.406 


± 


0.031 


4.056 


± 


0.280 


V497 Cep A 


6.89 


± 


0.46 


3.69 


± 


0.03 


4.290 


± 


0.021 


3.245 


± 


0.144 


V497 Cep B 


5.39 


± 


0.40 


2.92 


± 


0.03 


4.249 


± 


0.023 


2.877 


± 


0.152 



Table 4: Ages derived from fitting luminosities (^l) and radii (tn) of models for V906 Sco A and B to the 
observed values, for different combinations of hydrogen and heavy element abundances. While the mean 
of all values is 252 My, the mean value of t^A, ^lb and t^A is 256 My. 



X 


Z 


tMMy) 


tLB{My) 


tRA{My) 


twiMy) 


0.720 


0.0170 


256 


256 


253 


235 


0.729 


0.0156 


256 


259 


258 


240 


0.737 


0.0142 


257 


260 


253 


244 
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Prom the models of the component stars with solar 
composition {Xo = 0.70975, Zo = 0.016), we find the 
agreement time for the quantities La , Lb , Ra , and _Rb • 
While tLA and Ira arc about 3.5 My, Ilb and Irb arc 
about three and two times larger than 3.5 My, respec- 
tively. The reason for the small radius of V578 Mon 
A or large radius of V578 Mon B may be the accuracy 
of the data. However, if this is not the case then this 
system has some interesting features which arc beyond 
the standard models, for example, existence of rapidly 
rotating core for V578 Mon A (Yddiz 2003, 2005) or 
rotational mixing (Meynet & Maeder 2000). Within 
the standard approach, the present data implies that 
V578 Mon B is more evolved than V578 Mon A. The 
compatible radius for V578 Mon A based to its mass 
is about 6 Rq, which is 15 % larger than the observed 
value. Assuming the same effective temperature as 
given by Hensberge et al. (2000), the increase in lumi- 
nosity of V578 Mon A due to the increase in its radius 
is about 30 %. For this radius and luminosity, the age 
of the system is found as 6 My from the models of 
V578 Mon A and B with Z=0.01 and X=0.710. 

Hensberge et al. (2000) shows that the position 
of V578 Mon A in HRD corresponds to 11 Mq. This 
may be a manifestation of the rapidly rotating cores 
of early type stars. 

3.4 V392 Car 

The masses of the components of V392 Car are very 
close to each other: Ma = 1.90 ± 0.02 M©, Mb = 
1.853 ±0.02 Mq (Debernardi & North 2001). In Table 
5, the ages found from the luminosities and radii for 
different combinations of X and Z are listed. For such 
systems with similar components, small differences due 
to uncertainties in the observed quantities may lead 
very different results in model computations. How- 
ever, the observed values of the radii of the component 
stars are very accurate and provide an opportunity to 
determine the age of the system very accurately. The 
models with X = 0.6837 and Z = 0.0178 predict very 
close ages for V392 Car A and B (126 and 122 My, 
respectively). The adopted value for the age, 124 My, 
is the mean of these values. Debernardi & North 
(2001) reported very similar metallicity for V392 Car, 
Z=0.018, using evolutionary tracks of Schaller et al. 
(1992). 

3.5 V1229 Tau (HD 23642) 

The fundamental properties of components of V1229 
Tau are determined from light and radial velocity curves 
by Munari et al. (2004), Southworth et al. (2005) and 
Groenewegen et al. (2007). The results of these stud- 
ies are very similar. Among the eclipsing binaries in 
open clusters considered, this system has the minimum 
value of 

Vm = log(J?A/J?B)/log(MA/MB). 

It is not possible to fit models of V1229 Tau A and 
B to the observed properties simultaneously. There- 
fore, we find age of the system by considering V1229 
Tau A and B separately. 



The age derived from the calibration of V1229 Tau 

A is 193 My for {X, Z) = (0.6946, 0.01575) and 222 
My for {X,Z) = (0.7104, 0.014). However, the age 
derived from the calibration of V1229 Tau B is 730 
My for (X,Z) = (0. 7046,0. 01305). This value is about 
five times greater than the value given by Kharchcnko 
et al. In comparison to the age found for V1229 Tau B, 
the ages found for V1229 Tau A are close to the ages 
of Kharchcnko et al. Southworth et al. (2005) finds 
an age about 125-175 My, close to the results from the 
calibration of V1229 Tau A. The age of V1229 Tau is 
listed as 193 My in Table 2. 



3.6 V453 Cyg 

Recently, precise fundamental properties of V453 Cyg 
A and B have been presented by Southworth, Maxted 
& Smalley (2004b). However, there is a contradic- 
tion in their results. According to the radii and effec- 
tive temperatures of the components, the luminosities 
must be lower than the published values; in their Tar- 
ble 11, they give 

\og{LA/LQ^) = 4.69 and log(LB/is)=4.24 but we com- 
pute the luminosities from the radii and effective tem- 
peratures given in that table and find that log(LA / Lq ) — 
4.513 and log(LB/i'o)=4.056. The correction in log(L) 
is about 0.18; about 50% (see Table 3). 

From the calibration of model radii of V453 Cyg A 
and B for X=0.700, we find the age of the system as 
10.3 My and Z=0.0174. For X=0.7046, however, the 
age and Z are computed as 10.5 My and 0.0117, re- 
spectively. Whereeis the ages found with very different 
chemical compositions are very close, the model lumi- 
nosities with the latter values arc in perfect agreement 
with the (corrected) values found from observations. 
Southworth, Maxted & Smalley (2004b) give the age 
of the system as 10 My. The ages determined for V453 
Cyg are in good agreement with this age and the age 
of the cluster given by Kharchenko et al. (2005). Such 
a good agreement is based on the fact that V453 Cyg 
A is very close to the TAMS. 

Abundance analysis of V453 Cyg A and B is made 
by Pavlovski & Southworth (2009). They found that 
N, O and Mg abundances for V453 Cyg A are nearly 

solar and C, Si and Al are undcrabundant. Regarding 
the fact that O is the most abundant heavy element 
for normal stars, we estimate that Z ~ 0.01 from oxy- 
gen abundance (-0.08 dex) estimated by Pavlovski & 
Southworth (2009). This value is in good agreement 
with the values predicted from the models (Z=0.0117). 
Pavlovski & Southworth also re-derived the effective 
temperatures of V453 Cyg A and B: Tcha = 27, 900 ± 
400 K and T^ffB = 26, 200 ± 500 K. These values are 5 
and 3% greater than the values given by Southworth, 
Maxted & Smalley (2004b). 

Pavlovski & Southworth forecast the helium abun- 
dance as 0.086 and 0.123 dex, depending upon the 
method for derivation of microturbulence velocity. For 
these helium abundances, very similar ages are found: 
for 0.086 dex, Z=0.0113 and the mean age is 10 My; 
for 0.123 dex, Z=0.0163 and the mean age is 9 My. 
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Table 5: Ages derived from fitting the luminosities {ti,) and radii (^r) of models for V392 Car A and B 
to the observed values, for different combinations of hydrogen and heavy element abundances 

X Z ti^AjMy) tLB(My) tRAjMy) tnBjMy) 



0.7200 


0.0135 


260 


200 


273 


284 


0.7020 


0.0155 


200 


156 


189 


203 


0.6900 


0.0174 


220 


182 


140 


156 


0.6837 


0.0178 


157 


98 


126 


122 


0.6825 


0.0180 


163 


101 


119 


116 


0.6880 


0.0184 


246 


182 


112 


105 


0.6800 


0.0184 


165 


99 


111 


102 



3.7 V497 Cep 

The fundamental properties of this system were ob- 
tained by Yakut et al. (2003). Using their findings, 
the best agreement between the models and the ob- 
servational results from photometric and spectroscopic 
analysis is reached with X=0.744 and Z=0.0125. The 
derived age 23.0 My is the average of tLA = a ~ 22.3 
My, fLB=24.9 My, and tRB=22.5 My. Whereas the 
age wc find is close to the age for the cluster given in 
WEBDA, it is half of the value derived by Kharchenko 
et al. (2005). 

The accurate values given in Yakut et al. (2003) 
do not exactly obey L — AirR^aT^. It is not easy to 
find the source of this inconsistency. However, it may 
be due to round off. 

3.8 V381 Car (HD 92024) 

V381 Car A is a /3 Ccp for which three oscillation fre- 
quencies are observed with 1=2 and 4 (Frcyhammcr et 
al., 2006). Despite the large uncertainty in the funda- 
mental properties of its components (Preyhammer et 
al., 2005), it is a normal system with its slopes for 
the mass-radius and mass-luminosity relations given 
as 0.861 and 3.4894. From the fit of model radii and 
luminosities to the observed values, with X=0.738 and 
Z=0.0258, we find that tLA = t^A = 10.8 My and 
tRB=9.8 My. The age of the system is found as 10.5 
My, the average of these values. This age is in agree- 
ment with the age range (10-13 My) estimated by Prey- 
hammer et al. (2005). 

The chemical composition found from fitting inte- 
rior models of V381 Car A and B is very different from 
that of the other binaries in the sense that both of X 
and Z are greater than the solar values. This may be 
a result of the over-estimated masses of V381 Car A 
and B. 

3.9 V1034 Sco 

It is a normal system with its slopes for the mass- 
radius and mass-luminosity relations of 0.9876 and 
3.5728, respectively. From the fit of model radii and 
luminosities to the observed values, with X=0.6961 
and Z=0.02255, we find that tnA = tne = 5.64 My, 
tLA=5.72 My and iLB=5.05 My. The age of the sys- 
tem is computed as 5.51 My by talcing average of these 
values. Bouzid et al. (2005) confirm that the age of 



the system is about 5 My, in good agreement with the 
findings of the present study. 

3.10 DW Car 

At first sight, this system seems a normal system with 
slopes for the mass-radius and mass luminosity rela- 
tions of 0.9120 and 4.9858, respectively (Southworth 
& Clausen, 2007). According to the radii, the com- 
ponents of DW Car are unevolved MS stars and very 
close the zero-age MS (ZAMS). However, the luminos- 
ity derived from the observation implies that DW Car 
A is an evolved MS star. No simultaneous solution 
is available for the luminosities of DW Car A and B. 
From their radii, the age of the system is calculated as 
5.22 My. The non-rotating models giving this age are 
constructed with X=0.6934 and Z=0.01063. However, 
DW Car A and B are fast rotators; UeqA = 182 ± 3 
km s^^ and v^qB = 177 ± 3 km s^^. Therefore, mod- 
els rotating like a solid body are also constructed for 
them. The calibration of rotating models results in 
nearly the same age, 5.31 My. The hydrogen and metal 
abundances for the rotating models are X=0.6951 and 
Z=0.0078, very close to that of non-rotating models. 
Southworth & Clausen (2007) reported similar results; 
about 6 My, Z ~ 0.01. 

3.11 GV Car 

This system consists of two metallic-line stars and it 
shows apsidal motion. From the calibration of models 
of the components, with X=0.7046 and Z=0.0104, we 
estimate that tRA =324 My, tw = 372 My, tLA= 295 
My and tLB=302 My. The age of the system is found 
as 322 My. Southworth & Clausen (2006) give a very 
similar age for GV Car, 360 My. 

4 Mass and Time Dependence 
of the Mass— radius and the 
Mass— luminosity Relations 

The luminosity and radius of a non-rotating star are 
function of its mass, chemical composition (X and Z), 
and age: 

L = L{M,X,Z,t), R = R{M,X, Z,t). 

Changes in X and Z shift the evolutionary track of 
a model for a given mass in HRD. Despite the variar 
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Table 6: Luminosity and radius difference between ZAMS and TAMS for evolutionary tracks of V906 
Sco A. Alog(£) = (Ltams - -t'ZAMs)/-t'TAMS, A log(i?) = (-Rtams - -Rzams)/-Rtams 



star 


X 


Z 


Alog(L) 


Alog(i?) 


V906 Sco A 


0.702 


0.0210 


0.458 


0.439 




0.720 


0.0170 


0.456 


0.444 




0.729 


0.0156 


0.457 


0.448 




0.737 


0.0142 


0.456 


0.446 


GV Car A 


0.705 


0.0200 


0.431 


0.422 




0.720 


0.0100 


0.442 


0.423 



4.5 



log(MA/MB) 



3.5 



II 
HD 23642 (M4=2.193 Mq 
V(il8 Per Mi=2.332 Mq 
- GV Car (M4=2.r)l Mq 
V615 Per (Ma =4.075 Mq 
V497 Cep (Ma=6.89 Mq 
V453 Cyg (Ma=14.36 Mq 
V578 Mon Ma=14.54 Mq 
-V1034 Sco (Ma=16.84 Mq 


1 1 1 1 








I I 


1 I 1 I 



0.1 0.2 0.3 



0.4 0.5 0.6 0.7 0.8 

t/tMSA 



Figure 2: The mass-luminosity relation for the 
DEBs in open clusters is plotted with respect to 
the relative age of the primary components. 
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V618 Per Ma=2.332 Mq 
GV Car (Ma=2.51 Mq 
V615 Per (Ma =4.075 Mq 
V497 Cep (Ma=6.89 Mq 
V453 Cyg (Ma=14.36 Mq 
V578 Mon (Ma=14.54 Mq 
V1034 Sco Ma=16.84 Mq 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 O.E 

t/tMSA 



tion in positions of ZAMS and TAMS points according 
to chosen values of X and Z, ratio of TAMS to ZAMS 
luminosity is independent of X, Z and stellar mass: the 
logarithmic luminosity difference, Alog(L) = Ltams — 
i^zAMs)/i'TAMS , is Constant. In Table 6, A log(L) and 
A log(J?) for the models of V906 Sco A and GV Car A 
are given in the fourth and the fifth columns, respec- 
tively. For very different combinations of X and Z, 
the differences are between 0.42 and 0.46: the ZAMS 
luminosity (radius) is about 44 % less than the TAMS 
luminosity (radius), at least for the early- type stars. 
Consider a binary system whose component masses are 
very different, so that the massive component is near 
the TAMS and the low- mass is still near the ZAMS. 
For such a system, the luminosity ratio {La /Lb) is at 
a maximum. It was at a minimum when both of the 
components were near the ZAMS, and gradually in- 
creased with time. This is the case also for the ratio 
of the radii of the components. In this section, we dis- 
cuss variations of radii and luminosity ratios with time 
and then develop two methods for age determination 
of eclipsing binaries. 

As stated above, the mass-luminosity relation is 
not constant and mainly depends on the stellar mass 
and it increases with time. In Fig. 2, the luminosity- 
mass ratio 

{Im = log(LA/iB)/log(A/A/MB)) obtained from tfic 
models with values given in Table 2 is plotted with 
respect to the relative age in MS (i/tMSA)- Whereas 
Im increases with time, there is an inverse relation be- 
tween Im and Ma- Near the ZAMS, Im varies between 
3 and 4 and takes values between 4 and 5 near the 
TAMS. In principle, Im is very useful for predicting 
the ago of a system if its accuracy is good enough. 

The MS life time of the primary stars Imsa, given 
in the seventh column of Table 2, is mainly a func- 
tion of Ma, despite very different chemical composi- 
tion found for the DEBs. From a fit to the Ma-^msa 
relation we find that 



3.37 X 10" 
(M/M0)2-i; 



(4) 



4.1 



Age Determination from the Slope 
of the Mass— radius Relation 



Figure 3: The mass-radius relation for the DEBs 
in open clusters. 



The slope of the mass-radius relation 
Vm = log(JiA/i?B)/log(MA/MB) is also a function of 
the primary mass and time. In Fig. 3, rm is plotted 
with respect to t/tMSA- From the curves in Fig. 3, 
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Figure 4: The ages of binaries derived froni the 
ratio of radii (filled circle) are plotted with respect 
to the ages found from fitting models. For com- 
parison, ages given by Kharchenko et al. (2005) 
are also shown (circle). 



we determine the relative age of the binary systems 
by using the observed values of rm- Multiplication of 
relative age of the primary component of a binary sys- 
tem by tmsA gives us age of that binary. The ages 
(irm) of all the binaries computed with this method 
arc listed in the eighth column of Tabic 2. In Fig. 4, 
these ages (hlled circles) arc plotted with respect to the 
ages from fitting interior models of component stars to 
the observed accurate values (the fifth column in Table 
2). For comparison, the ages found by Kharchenko et 
al. (2005) arc also shown (circle). There is in general 
a good agreement between the ages, except for V6I8 
Per and V578 Mon. The observed value of Vm is very 
high for V618 Per and low for V578 Mon, in compar- 
ison with the model values. For V578 Mon, Garcia et 
al. (2010) report that they find a different value for 
the secondary radius. The difference may significantly 
influence the observed value of rm- 

The age of V1229 Tau is computed as 124 My by 
using Vm from Munari et al. (2004). rm from South- 
worth et al. (2005) gives age of the system as 97 My. 

The precision of the ages found from the mass- 
radius curves depends on how sensitive the relative age 
is to chemical composition. In Fig. 5, r^s derived from 
models of V618 Per A and B with different chemical 
compositions are plotted with respect to the relative 
age t/tMSA. For X=0.705, there is a perfect agreement 
between r^s of models with Z=0.01 (thin solid line) 
and Z=0.02 (dotted line). We shall also compare Vm 
values of models with the same Z but different X. Vm 
of models with X=0.735 and Z=0.01 (thick sohd fine) 
is slightly less than that of the models with X=0.705 
and Z=0.01. The agreement between the three curves 
of r-m is good enough to suggest that relation between 



Figure 5: The mass-radius relation for models of 
V618 Per with different X and Z. 

Vm and relative age relation is almost independent of 
chemical composition. For the precision of the ages, 
the sensitivity of iMSA to chemical composition is also 
important. As given in Eq. 4, ^msa is a very strong 
function of stellar mass. The secondary effect of Z on 
tmsA behaves as tmsA oc (Z/0.015)° ''^ and is neglected 
in the present study. 

4.2 Simple Expressions for the r^- 
age Relation 

The time dependence of rm of binary systems with very 
different masses and chemical compositions is very sim- 
ilar and therefore can be used to estimate the age of a 
binary system. To do this, the time variation of rm is 
fitted with an equation of the form 

Vm = a(t/tMSA)" + b. (5) 

According to Ma of the binaries, the curves of Vm 
in Fig. 3 separate essentially into two groups. The r-m 
values of the binaries with Ma <3.4 M© take place in 
the lower part, and those with Ma > 3.4 M© appear 
in the upper part of Fig. 3. Near the ZAMS of the 
primary component, rm is about 0.4 for Ma < 3.4 Mq 
and 0.6 for Ma greater than 3.4 Mq. It increases with 
time and reaches a value between 1.6 and 2.4. It seems 
that rm is a better indicator of binary system age than 
Im- For example, age of a system with = 1 is about 
half of the TAMS age of the primary component. For 
more accurate ages, two fits for the two mass ranges 
are made to the curves in Fig.3. 

4.2.1 DEBs with Relatively Low- mass Pri- 
mary Component (Ma < S.AMq) 

In order to increase the precision of Eq. (5), we derive 
two expressions for low and high values of Vm- 
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For the three binary systems with Ma < 3.4Mq, 
namely, V618 Per, V1229 Tau (HD23642) and GV Car, 
three separate fits are made to the curves in Fig. 3 by 
using Eq. 5. Then, the mean values of the parameters 
(a, b and n) given in Eq. 5 yield 

— ^ = 0.941 (r^ - 0.39)° ®^* for < 1.0 (6) 

IMSA 

and 

= 0.77{rm - 0.588)°-**' for Vm > 1.0. (7) 

IMSA 

The ages of V1229 Tau and GV Car computed from 
these expressions (trmf) are given in the ninth column 
of Table 2. They are in good agreement with the ages 
found by other methods. 

4.2.2 DEBs with Relatively High-mass Pri- 
mary Component (Ma > 3AMq) 

Using Eq. 5, two fits are applied to the curves of V615 
Per, V497 Cep and V453 Cyg given in Fig. 3 according 
to the value of r^- We compute the average values of 

the parameters a, b and n. From these values, the 
expressions for age as a function of rm are obtained: 

— ^ = 1.011(r,„ - 0.561)°-^^^ for rm < 1-0 (8) 

IMSA 

and 

— ^ = 0.613{rm - 0.852)°-^^° for > 1.0. (9) 

tmsA 

For V615 Per, rm = 0.747 and Eq. (8) gives its age as 
27.1 My. These two expressions are used to determine 
age of the binary systems whose Ma greater than 3.4 
Mo (see Table 2). 

The fitting formula given in Eqs. 6-8 arc not valid 
for V381 Car because of the large mass difference be- 
tween its components. They are also not applicable to 
binaries with very similar components. For V392 Car, 
for example, the mass (2.47%) and radius (1.54%) dif- 
ferences between its components are so low that the 
fitting formula gives an age two times greater than the 
ages found by other methods. The situation for DW 
Car is better than V392 Car. The mass (6.26%) and 
radius (5.73%) differences between the components of 
DW Car are big enough to predict a reasonable value 
for age. For V906 Sco, however, there is a big differ- 
ence between the radii of the components which gives 
so great value for rm that we deduce that the age is 
about trmf = tMSA ~ 253My. 

The rm curve of V1034 Sco (x in Fig. 3), which 
has the highest primary mass (Ma= I6.84M0) of the 
stars studied, is slightly different from the rest. It is 
in the upper part during the early MS phase as that 
of the other binaries with Ma > 3.4 Mq, and joins 
later on the curves for binaries with Ma <3.4 Mq. 
Therefore, the expressions given in Eqs. 6-7 shoul be 
used rather than Eqs. 8-9 for Ma > 15 Mq when the 
relative age is greater than 0.5. 

The new expressions proposed above were addi- 
tionally tested for selected eclipsing binaries. The age 
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Figure 6: The mass-radius relation for the DEBs 
given in Andersen (1991) with respect to the rela- 
tive age of the primary components. The ages are 
taken from Pols et al. (1997). The solid line is for 
the curve fitted to this data by using Eq. 5. For 
comparison, of V497 Cep from models (dotted 
line) is also plotted. 

of HD 42401 is given by Williams (2009) as 25 My 
and we find the same age ( 25.6 My) from the ratio of 
radii. For the age of V539 Ara Clausen et al. (1996) 
and Lastennet and Valls-Gabaud (2002) give 45 and 
39 My, respectively. The ratio of radii gives the age as 
40.0 My. 

4.3 Age Determination from the Slope 
of the Mciss-luminosity Relation 

We also derive an expression to determine age from 
the ratio of luminosities: 

t = 0.7M{lm - 4.535Ma" '''*)10.('°'^*'*-°-''' Imsa- 

(10) 

This expression gives acceptable ages for only the five 
binary systems (V453 Cyg, V497 Cep, V381 Car, V1034 
Sco and GV Car). The reasons might be that the ra- 
tio of luminosities of component stars derived from the 
light curve of a binary system may not be as accurate 
as the ratio of radii. The dependence of Im on Ma is 
also as important as its time dependence (see Fig. 2). 

4.4 Does an r^— relative Age Rela- 
tion Implicitly Exist in the Lit- 
erature? 

The time dependence of rm is not considered in the 
previous studies in the literature. However, this de- 
pendence also exists in Pols et al. (1997), for example, 
though the authors are not aware of it. In Fig. 6, 
rm is plotted with respect to the ages they find (in 
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units of MS lifetime of the primary components) for 
the well known binaries (Andersen 1991). In order to 
avoid complications, the binaries with Ma > 1.5 and 

significantly different components (difference between 
masses of the components is greater than 10%) are 
considered. The mass range is 1.56-9.245 M©. The 
solid line shows the fitted curve for this data using Eq. 
5. The dotted line represents the rm of V497 Cep, al- 
ready plotted in Fig. 3. The agreement between the 
two curves is amazing. 

5 The Uncertainty in the Ages 
of the Binaries 

The ages of stars we have found involve many uncer- 
tainties, based on the model ingredients (e.g., opac- 
ity, nuclear reaction rates, equation of state) and the 
observational quantities. It is very difficult to com- 
pute the uncertainty duo to the former and therefore 
we take care in calculating the uncertainties resulting 
from the latter in our error analysis. The uncertainty 
in the ages (At) derived from the luminosities and radii 
of component stars is the time interval in which evo- 
lutionary tracks of both components stay within the 
corresponding error boxes in luminosity-radius (L-R) 
diagram. The values of At for the binaries are given 
in the third column of Table 7. For a binary system 
whose age is derived from its primary component, the 
uncertainty is also computed from that star. Within 
the context of our approach, the uncertainties are rear 
sonably small. 

We also compute uncertainties for the ages (Atrm) 
derived from the ratio of the radii {tab ~ Ra/Rb)- 
The uncertainty in the mass ratio appearing in the de- 
nominator of Trm is usually very small for well-measured 
binary systems. The values of ArAB are available in 
the literature and given in the fifth column of Table 
7. The error in age due to uncertainty in the ratio of 
the radii is given in the seventh column of Table 7. 
The uncertainty in the relative age {Atrei = At/tMSA) 
derived from the ratio of the radii (rAB = Ra/Rb) is 
computed using 



Atrel = 



ArAB 



(drAB/dtrel) ' 



(11) 



The derivative of the ratio of radii with respect to the 
relative age in Eq. 11 is obtained from the models 
around the relative age of each binary system. As the 
masses of component stars are close to each other the 
derivative is negligibly small and therefore the uncer- 
tainty in age is very large. This is the case for the 
binary systems whose components are very similar: 
V906 Sco, V392 Car and DW Car. Except V615, the 
rest of the binary systems have very small uncertain- 
ties for their ages. For V615 Per, the reason for the 
large uncertainty is the large value of ArAB. 

6 Conclusions 

Our current understanding of the evolution of our far 
and near universe depends on accurate age estimates 



for astrophysical systems. The stars are the main con- 
stitutes of these systems, whose ages are essentially de- 
termined from comparison of their observed and model 

properties. Rotating (for DW Car A and B) and iiou- 
rotating models for the components of the eclipsing 
binaries, members of open clusters, are constructed to 
derive the ages and chemical compositions of the bi- 
naries and hence the clusters. The largest difference 
between the ages we found for the eclipsing binaries 
and the ages given by Kharchenko et al, (2005) for 
the clusters occurs for the V497 Cep binary system in 
NGC 7160, where the difference is about 100 %. For 
the remaining systems the difference is less than 40%. 

The binary systems such as V906 Sco in NGC6475 
and V453 Cyg in NGC 6871 are the most suitable 
systems for age determination: their primary compo- 
nents are about to complete their MS life time and the 
secondary components are near or not far from the 
ZAMS. For such systems, despite equivalent solutions 
with very different chemical compositions, the ages are 
very close. For the three different solutions with metal- 
licity 0.014-0.017, the age of V906 Sco is in the range 
252-254 My. 

The slope of the mass-luminosity relation for a 
binary system depends on mass of its primary compo- 
nent Ma, if Ma is less than about 10 M©. The value 
of 

= log (L A /ie)/ log (Ma /Mb) near the ZAMS is 4 
for Ma=2.5 Mq, 3.5 for Ma=6.89 Mq and it drops 
to a minimum value of 3 for Ma=10 M0. For binary 
systems with Ma greater than 10. M©, = 3 near the 
ZAMS and reaches to 4.3 near the TAMS age of their 
primary component. 

The slope of the mass-radius relation changes 
with time more meaningfully than the slope of the 
mass-luminosity relation, so is more suitable for esti- 
mating the age of the system once it is calibrated. If 
rm=l for a system, its age is about half of the MS 
lifetime of its primary (massive) component. We have 
derived simple expressions to estimate the age of a bi- 
nary system from the slope of mass-radius and mass- 
luminosity relations. In particular, the expressions for 
the slope of mass-radius relation give very reasonable 
results (irm). Eqs.(6)-(9) can be used to predict age of 
any binary system if the radii and masses of its com- 
ponents are known. The new expressions are tested 
for the eclipsing binaries HD 42401 and V539 Ara and 
give excellent results for their ages. The advantage of 
these expressions for the relative age {t/tyisA) is that 
they are independent of chemical composition. 
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